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ABSTRACT 

We study the acceleration of heavy nuclei at SNR shocks taking into account the 
process of ionization. In the interstellar medium atoms heavier then hydrogen which 
start the diffusive shock acceleration (DSA) are never fully ionized at the moment 
of injection. We will show that electrons in the atomic shells are stripped during 
the acceleration process, when the atoms already move relativistically. For typical 
environment around SNRs the dominant ionization process is the photo-ionization due 
to the background Galactic radiation. The ionization has two interesting consequences. 
First, because the total photo-ionization time is comparable to the beginning of the 
Sedov- Taylor phase, the maximum energy which ions can achieve is smaller than the 
standard result of the DSA, which predict -E m ax oc Z^. As a consequence the structure 
of the CR spectrum in the knee region can be affected. The second consequence is that 
electrons are stripped from atoms when they already move relativistically hence they 
can start the DSA without any pre-acceleration mechanism. We use the linear quasi- 
stationary approach to compute the spectrum of ions and electrons accelerated after 
being stripped. We show that the number of these secondary electrons is enough to 
account for the synchrotron radiation observed from young SNRs, if the amplification 
of the magnetic field occurs. 

Key words: shock acceleration - ionization - non thermal emission - supernova 
remnant 



1 INTRODUCTION 

The bulk of Galactic cosmic rays (CRs) is largely thought 
to be accelerated at the shock waves associated with super- 
nova remnants (SNRs) through the mechanism of diffusive 
shock acceleration (DSA). A key feature of this mechanism 
is that the acceleration rate is proportional to the particle 
charge. This property is especially appealing if one tries to 
explain the structure of the knee in the CR spectrum: if we 
assume that the maximum energy of accelerated particles 
scales with the charge of the particles involved, a knee arises 
naturally as a superposition of spectra of che micals with dif- 
ferent nuclear charges Zjsie |Horandej|2003h . This result is 
based on the assumption that nuclei are completely ionized 
during the acceleration. On the other hand when ions are 
injected into the acceleration process they are unlikely to be 
fully stripped, especially if of high nuclear charge. 

The atoms relevant for the injection are those present 
in the circumstellar medium where the forward shock prop- 
agates. The temperature of this plasma varies from 10 4 K 
if the SNR expands into the regular ISM, up to 10 6 K if 
the expansion occurs into the bubble created by the pro- 
genitor's wind. If T ~ 10 4 K even hydrogen is not fully ion- 
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ized, as demonstrated by the presence of Balmer lines associ- 
ated with shocks in some young SN Rs (|Chevalier et aljl98d : 
ISollerman at al.ll2003l : lHendl2009l ). For T ~ 1 6 only atoms 
up to Zn = 5 can be completely ionized (|Porquet et al.l 
l200ll ). 

The typical assumption made in the literature is that 
atoms lose all electrons in the atomic orbitals soon after the 
beginning of the acceleration process, namely that the ion- 
ization time needed to strip all the electrons is much smaller 
than the accel eration time. In spite of this assumption in 
iMorlind (2009) we showed that, for a typical SNR shock, 
the ionization time is comparable with the acceleration time, 
hence electrons are stripped when ions already move rela- 
tivistically. This fact has two important consequences: 1) the 
maximum energy of ions can be reduced with respect to the 
standard prediction of DSA and 2) the ejected electrons can 
easily start the acceleration process because they already 
move relativistically. 

The possibility that ionization can provide a source of 
relativistic electrons is especially relevant because the ques- 
tion of how electrons are injected into the DSA is still an 
unsolved issue. DSA applies only for particles with a Lar- 
mor radius larger than the typical shock thickness, which is 
of the order of Larmor radius of the shocked downstream 
thermal ions. The injection condition can be easily fulfilled 
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for supra-thermal protons which reside in the highest energy 
tail of the Maxwellian distribution. The injection of heavier 
ions is, in principle, even easier than that of protons in that 
their Larmor radius is larger assuming they have the same 
proton temperature (which could be not the case). On the 
other hand the same argument tells that electrons cannot be 
injected from the thermal bath because, even if we assume 
equilibration between electrons and protons, the Larmor ra- 
dius of electrons is a factor (m e /m p ) 1//2 smaller than that 
of protons. Only electrons which are already relativistic can 
cross the shock and start the DSA. 

Most proposed solutions to the electron injection prob- 
lem involve some kind of pre acceleration mechanism driven 
by plasma instabilities able to accelerate therm al elec- 
trons up to mildly relativistic energy. For instance iGaleevi 
(| 1984T) showed that magnetosonic turbulence excited by 
a beam of reflected ions ahead of the shock can acceler- 
ate electrons along the magnetic field lines thanks to their 
Cerenkov interaction with excited waves. Other studies pre- 
dict that electrons can be effectively preaccelera ted by self- 
gener ated whistler waves ()Levinsonlll992l . ll996l ). More re- 
cently |Xmano^&jIoshin3 l|2009h showed that efficient elec- 
tron acceleration can happen in perpendicular shocks due to 
"shock surfing" of electrons on electrostatic waves excited by 
Buneman instability at the leading edge of the shock foot. 
A diff erent mechanism proposed bv lRiquelme fc Spitkovskvl 
i|2010l) . can take place in quasi-perpendicular shocks and is 
driven by the presence of oblique whistler waves excited by 
the returning ions in the shock foot. The study of these elec- 
tromagnetic pre-acceleration mechanisms can be only per- 
formed using numerical techniques which showed that the 
fraction of injected electrons strongly depend on the values 
of initial conditions (like magnetic field strength and orien- 
tation) and are difficult to apply to realistic cases. 

In this paper we investigate in detail the process of ion- 
ization applied to the acceleration of heavy ions, presenting 
the full steady-state solution for shock acceleration in the 
test particle approach. 

The paper is organized as follows: in Ej2] we compare 
the acceleration time of ions with the ionization time due 
both to photo-ionization and to Coulomb scattering show- 
ing that the former dominates on the latter. In ij3]we use the 
linear acceleration theory to compute the maximum energy 
achieved by different chemical specie when the ionization 
process is taken into account. In [|4]we compute the distri- 
bution function of both ions and electrons in the framework 
of linear acceleration theory including the term due to ion- 
ization. We conclude in i|5] 



2 IONIZATION VS. ACCELERATION TIME 

In this section we show that the ionization of different chem- 
ical species occurs during the acceleration process on a 
timescale which is comparable with the acceleration time 
needed to achieve relativistic energies. We also show that 
the ionization due to Coulomb collision is generally negligi- 
ble compared with the photo-ionization due to the interstel- 
lar radiation field (ISRF) . 

Let us consider atoms of a single chemical specie N with 
nuclear charge Zn and mass mjv = Am p , which start the 
DSA with initial charge Z < Zn and momentum Pmj. We 



want to compute which is the ionization time needed to lose 
one electron, changing the net charge from Z to Z + 1 and 
which is the momentum p that ions reach when ionization 
occurs. 

For simplicity we compute the acceleration time in 
the framework of linear shock acceleration theory for plane 
shock geometry: i.e. we assume that during the ionization 
time needed to strip one single electron, the shock structure 
does not change. If a particle with momentum p diffuses 
with a diffusion coefficient D(p), the well know expression 
for the acceleration time is: 

/ U! — U 2 \ U! U2 P 

" ^liij N ' 

where u is the plasma speed in the shock rest frame, and the 
subscript 1 (2) refers to the upstream (downstream) quanti- 
ties (note that « s hock = ^i). The downstream plasma speed 
is related to the upstream one through the compression fac- 
tor r, i.e. U2 = ui/r. We limit our considerations to strong 
shocks, which have compression factor r — 4, and we assume 
Bohm diffusion coefficient, i.e. Db = rL/3c/3, where /3c is the 
particle speed and tl = pc/ZeB is the Larmor radius. The 
turbulent magnetic field responsible for the particle diffu- 
sion is assumed to be compressed downstream according to 
B2 = rBi. Even if this relation applies only for the magnetic 
component parallel to the shock plane, such assumption does 
not affect strongly our results. It is useful to define the in- 
stantaneous acceleration time as t aC c = dr /(dp/p). Using all 
previous assumptions t acc can be expressed as follows: 

t acc (p) = 0.85 us 2 (£) yr . (2) 

Here the upstream magnetic field is expressed in fiG and the 
shock speed is ui = u$ I0 8 cm/s. Notice that, when B and 
u are assumed constant, Tacc reduces to tacc for pi n j <S p 
and D tx p. To compute the energy reached by particles 
when the ionization event occurs, Eq. ((2)1 has to be com- 
pared with the ionization time scale. Ionization can occur 
either via Coulomb collisions with thermal particles or via 
photo-ionization with background photons. Whether the for- 
mer process dominates on the latter depends on the ISM 
number density compared to the ionizing photon density. In 
the following we show that photo-ionization dominates when 
SNRs are young and expand in a typical ISM. 

As soon as the ions start the DSA, the photo-ionization 
can occur only when the energy of background photons, e', as 
seen in the ion rest frame, is larger than the ionization energy 
I. Atoms moving relativistically with a Lorentz factor 7 see 
a distribution of photons peaked in the forward direction 
of motion, with a mean photon energy e' = 'ye. The photo- 
ionization cross section can be estimated using the simplest 
approximation for the _K"-shell cross section of hydrogen- like 
atoms with effective nuclear charge Z, i.e. (jHeitlerll 19541 ): 

Oph(e') = 64a~ 3 a T Z- 2 (l N ,z/e') 7/2 (3) 

where ot is the Thompson cross section and a is the fine 
structure constant and In,z is the ionization energy thresh- 
old for the ground state of the chemical specie N with Zn — Z 
electrons. The nu merica l value s of In,z can be found in the 
literature (see e.g. lAllenl (|l973l ) for elements up to Zn = 30) . 
To get the full photo-ionization time we need to integrate 
over the total photon energy spectrum, i.e.: 
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^(7) = j 'de^^-ca^e), (4) 

where dn p h /de is the photon spectrum as seen in the plasma 
rest frame. Because the photo-ionization cross section de- 
creases rapidly with increasing photon energy, for a fixed 
ion speed the relevant ionizing photons are only those with 
energy close to the threshold, i.e. e ~ In,zH, measured in 
the plasma frame. The corresponding numerical value is: 

r ph ( 7 ) ~ 0.01 Z 2 (n ph (Jiv,z/7)/cm- 3 ) ~* yr . (5) 

We assume that the maximum possible acceleration time 
is equal to the Sedov- Taylor time, tsr, corresponding to 
the end of the free expansion phase. Comparing tsr with 
Eq. ([5]), we see that the photons which can be relevant 
for ionization are only those with a number density n pn > 
10 _5 Z 2 (tsr/10 3 yr) _1 cm -3 . For this reason we can neglect 
the high energy radiation coming from the remnant itself 
because the typical number density of the x-ray photons is 
less than 10 -7 ph/cm -3 . 

An important consequence of the photo-ionization is 
that electrons are stripped from the atoms with a kinetic 
energy <C m e c 2 as measured in the atom rest frame. Hence, 
if the parent atoms move relativistically in the plasma rest 
frame, electrons move with the same Lorentz factor and the 
same direction of the parent atoms. This property will be 
used in 33]to simplify the calculation of electrons spectrum. 

Now we consider the ionization of accelerated ions due 
to Coulomb collisions with protons and electrons of the ther- 
mal plasma. We neglect the thermal energy of those particles 
and we consider the process in the ion's rest frame. The role 
of collisions is mainly relevant during the very initial stage 
of acceleration, namely when the kinetic energy of incident 
particle, Ukin, is close to the ionization energy, /, as seen in 
the ion's rest frame. This is a consequence of the shape of 
ionization cross section which has a peak value for 15km ~ I 
and rapidly decreases for larger energies. A good upper limit 
for the peak value of the cross section is pro vided by th e 
classical (non relativistic) Thomson approach (|Allenlll973T ). 
which gives a maximum for Skin = 21 equal to: 

max 2 A7 j — 2 / c\ 

Cecil = ira N e I Ryd . (6) 

where ao is the Bohr radius, Inyd is the ionization energy 
expressed in Rydberg unity and N e is the number of elec- 
trons in the atomic orbital considered. For larger energies 
the Thomson cross section decreases like E^ while the 
asymptotic behavior predicted by the quantum mechanics is 
a ~ log(_Ekin)/£ : kin. This implies that the collisional ioniza- 
tion time has a minimum for a — cr™^ 1 and than increases. 
In order to get this minimum value one must average the 
contribution of collisions in the upstream and downstream 
region. Assuming equal density for thermal proto ns and elec- 
trons n e = %= f>i, the result is (|Morl ino 2009): 

rSir w [ca^u m(l + r)]- 1 = 0.0024 7 Ryd (^b) yr -(7) 

Equation J7J is valid in the non relativistic regime. In this 
regime the acceleration time increases only like t acc ~ p ~ 
E^ 2 . As a consequence the collisional ionization can occur 
only if the acceleration time needed to have -Ekin = 21 is 
larger than t™}". Let us call 7* the Lorentz factor which 
corresponds to this condition, namely 7* = 1 + 2I/m e c 2 . 



Comparing Eq. with Eq. ((2J, the condition t a cc(7*) = 
T coU provides a lower limit for the upstream density, i.e.: 

m.min = 277Ryd-B M GMg (Z/Zn) cm -3 . (8) 

Hence for m < ni.min ions are accelerated up to relativistic 
energies before being ionized by collisions. Notice that we 
neglected collisions with atoms heavier than hydrogen: such 
contribution could reduce the ionization time by less than a 
factor 2, hence are main conclusions remain valid. 

Now, in order to understand the relevance of collisions 
in the relativistic regime we can use the asy mptotic Beth e 
cross section valid in the limit -Ekm S> m e c 2 , (|Bethelll930n . 
The relativistic formulation of this cross section for the total 
ionization is often written as (|Kim at aUl200Ch : 

<TRBehte(7) = [M 2 (ln( 7 2 /3 2 ) - + C R ] . (9) 

The two constant M 2 and Cr are related to the atomic form 
factors of the target and are independent of the incident 
particle energy. Their exact numerical values are very diffi- 
cult to compute with some noticeable exception like H-like 
atoms. In principle they can also be inferred from photo- 
ionization experiments but in practice they are known only 
for a bunch of targets. Nevertheless from th eory we know 
that M 2 ~ NJ I Rvd , while C R ~ 10M 2 (|lnokutil Il97~il : 
iKim at alJl200Ch . Using these estimates in the ultra rela- 
tivistic limit we have CTRBethe/e>"™ii x ~ 10 -3 . This result im- 
plies that Coulomb collisions are much less important in the 
relativistic regime with respect to the non relativistic one. 

In order to compare the photo-ionization with the col- 
lisional ionization times, in Fig. Q] they have been plotted 
together with the acceleration time for three different H-like 
ions: He + , C 5+ and Fe 25+ . The characteristic time are plot- 
ted as functions of the ion's Lorentz factor. The acceleration 
time iacc, plotted with solid lines, is shown for two different 
choices of parameters: ui = 1000 km/s ; Bi = 3/iG (upper 
line) and for ui = 10 4 km/s, Bi = 20/xG (lower line). The 
photo-ionization time r p h (dashed lines) is computed accord- 
ing to Eq. Q using the cosmic mi crowave background plu s 
the Galactic ISRF as reported by iPorter fc Strong! (|2005T) . 
The ISRF strongly changes going from the center towards 
the periphery of the Galaxy, hence, the photo-ionization de- 
pends on the SNR location. In order to evaluate this vari- 
ation we show T p h using the ISRF in the Galactic Center 
(lower dashed line) and the one at 12 kpc far away from 
the Galactic Center in the Galactic plane (upper dashed 
lines). Finally with dot-dashed lines we plot the collisional 
ionization time for the upstream density ni = 1cm -3 : the 
thin horizontal line corresponds to the non relativistic lower 
limit expressed in Eq. T" C o!r> while the thick line is cal- 
culated using the relativistic Bethe cross section for H-like 
atoms. In this case the values of the parameters used in the 
Bethe cross section are M 2 = 0.30/Z% and Cr = 4.30/Z 4 ^ 
l|Bethelll930h . Fig.Q]clearly shows that the collisional ioniza- 
tion can be neglected for the case of young SNRs expanding 
into a medium with a number density less than a few cm -3 . 
The photo-ionization is the dominant process even for SNR 
located far away from the Galactic Center. This considera- 
tion is strengthened in the case of core-collapse SNRs, which 
typically expand into the bubble created by the progeni- 
tor wind, whose density is usually assumed ~ 10 _2 cm -3 . 
The only phase where the collisional ionization could domi- 
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nate over photo-ionization is the very initial stage, when the 
shock is expanding into the progenitor's wind, whose density 
is much higher. This phase typically lasts few years and in 
this work it will be neglected. 

In Fig.[T]the value of 7 where i acc cross r p h, identifies the 
Lorentz factor of the ejected electrons they are for. Looking 
at the upper plot it is clear that even electrons coming from 
the ionization of He + can, in principle, start the PS A, having 
a typi cal Lorentz factor > 10. In fact, as discussed in lMorlinq 
(2009) the minimum Lorentz factor required for electrons to 
be injected is -yi n j ~ 3 — 30 for typical parameters of SNR 
shocks. 



3 MAXIMUM ENERGY OF IONS 

The process of ionization can affect the maximum energy 
that nuclei achieve during the acceleration, especially those 
with large nuclear charge. It is worth stressing that the 
knowledge of the maximum energy is intimately connected 
with two important aspects of the CR spectrum: the in- 
terpretation of the knee structure and the transition region 
from Galactic to extragalactic component. The knee is com- 
monly interpreted as due to the superposition of the spec- 
tra of all chemicals with different cutoff energies. Using the 
flux of different components measured at low energies it has 
been shown that the knee structure is well reproduced if one 
assumes that the maximum energy of eac h specie -E max ,i v 
is proportional to the nuclear charge Zm (Horandel 2003). 
Nevertheless the superposition of subsequent cutoff seems 
to be confirmed by the measurements of the spectrum of 
single components in the knee region: data presented by 
the KASCADE experiments show that the maximum en- 
ergy of He is ~ 2 times larger than that of the protons 
jAntoni et al J 120051 ) . From the theoretical point of view the 
relation S ma x,AT oc Zm is clearly predicted by DSA if one as- 
sume that the diffusion coefficient is rigidity-dependent. On 
the other hand this is correct only provided that during the 
acceleration process nuclei are completely ionized otherwise 
the maximum energy is proportional to the effective charge 
of the ions rather that to their nuclear charge. 

Now, the maximum energy is thought to be achieved 
at the beginning of the Sedov- Taylor phase (tsr) 
jBlasi. Amato fc Caprioll 120071 ). For later times, t > tsr, 
the shock speed decreases faster than the diffusion veloc- 
ity, hence particles at the maximum energy can escape from 
the accelerator and the maximum energy cannot increase 
further. It is worth noting that the process of escaping 
is not completely understood in that it strongly depends 
on the turbulence produced by the same escaping parti- 
cles and on the damping mechanisms of the turbulence it- 
self (|Caprioli et al.| [2010). Anyway here we assume that the 
maximum energy, J5 max , is achieved at t — tsr- If the total 
ionization time is comparable, or even larger than tsr, we 
do expect that E max < E^a*, where we call -Emax the max- 
imum energy achieved by ions which are completely ionized 
since the beginning of acceleration. 

A consistent treatment of the ionization effects would 
require the use of time-dependent calculation. However our 
aim is to get a first order approximation of the effect of 
ionization. This can be achieved using the quasi-stationary 
version of the linear acceleration theory. 




Log(y) 

Figure 1. Comparison between acceleration time (solid lines), 
photo-ionization time (dashed) and collisional ionization time 
(dot-dashed) for three hydrogen-like ions: He+, C 5 + and Fe 25+ 
(from top to bottom). Acceleration and photo-ionization times 
are shown with two curves representing two different set of pa- 
rameters, as explained in the text. The horizontal thin dot-dashed 
line represent the lower limit for the collisional ionization time, 
Eq. JBJ, while the thick dot-dashed line is computed using the 
Bethe cross section with a plasma density ni = 1cm -3 . 



In order to compute E max the first piece of information 
we need is the initial ionization level of ions. In the case 
of volatile elements, which exist mainly in the gas-phase, 
the level of ionization is easy to estimate being only a func- 
tion of the plasma temperature: ions are ionized up to a 
level such that the ionization energy is of the same order 
of the kinetic energy of thermal electrons as seen in the 
ion's rest frame. On the other hand refractory elements are 
largely locked into solid dust grains. This occurs both in 
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the regular ISM as well as in stellar winds. As claimed by 
Ell ison, Drurv fc Meveil (|l997f) there are evidences that ele- 
ments condensed in dust grains are efficiently injected into 
the DSA thanks to the fact that dust grains can be eas- 
ily accelerated up to energies where the grains start to be 
sputtered. A non negligible fraction of the atoms expelled 
by sputtering are energetic enough to start the DSA. In 
this scenario the ionization level of ejected atoms is very 
difficult to predict in that depends on the structure of the 
grains. However a reasonable assumption could be that ions 
preserve all the electrons in the inner shells which are not 
shared in the orbitals of the crystal structure of the grains. 
Moreover, as noted bv lEllison. Drurv fc Meverl l|l997h . if the 
ejected ion is highly ionized its electron affinity is strong and 
electron-exchange with the atoms of thermal plasma could 
reduce the level of ionization up to the equilibrium one. For 
these reasons, in the following calculation, we neglect the 
complication arising from the dust sputtering process and 
we assume that the level of ionization is only determined by 
the plasma temperature. 

Let us consider a single ion injected with momentum 
po and total charge Z\. The ion undergoes acceleration at 
a constant rate oc Z\ during a time equal to the ionization 
time needed to lose one electron, r p h,i, when it achieves 
the momentum pi. Notice that we include only the photo- 
ionization because, as we saw in the previous paragraph, 
for typical density of the external medium around the SNR 
(no ~ 1cm -3 ), the Coulomb collisions can be neglected. In 
order to compute simultaneously r p h,i and p\ we have to 
equate the acceleration time with the ionization time: 

Tph,l(pi/m N c) = Tacc(P(hPl) = iacc(pi) — tacc(po) • (10) 

The last equality holds because we use Bohm diffusion and 
because we assume that the shock velocity and the magnetic 
field both remain constant during r p h,i. Moreover we saw 
that the photo-ionization occurs when ions already move 
relativistically, hence using Eq. @ we set f3 = 1. Eq. (| 10p 
gives 

pi po Z x B x u\ 

= 1- -, y 7 Tph.i (pi m N c) , (11 

m v c m p c 1.7 Zn 

where the subscript 1 label the quantities during the time 
Tph.i. In Eq. (|11[) B and u are expressed in units of nG and 
ug while T p h,i is expressed in yr. Once the background pho- 
ton distribution is known Eq. (|ll|l can be solved numerically 
in order to get pi and r p h,i. After the first ionization the ac- 
celeration proceeds at a rate proportional to Z2 = Z\ + 1 
during a time T p h,2, which is the time needed to lose the 
second electron. Applying Eq. (|11[) repeatedly for all sub- 
sequent ionization steps, we get the momentum when the 
ionization is complete. If the total ionization time, t^, is 
smaller than the Sedov time, in order to get the maximum 
momentum we need to add the further acceleration during 
the time (tsr — T^n)- The final expression for the maximum 
momentum can be written as follows: 

Zn-Zq 2 M 2 

= > TTt h0 (tsr — T ion J > , (12) 

m„c t-^t 1.7 Zn ^— ' 1.7 M 

k = l i=l 

where we neglected the contribution of po- The last term has 
been written as a sum over M time-steps in order to handle 
the case where magnetic field and shock speed change with 
time. 



Now we can quantify the effect of ionization in the de- 
termination of p max . The simplest approximation we can do 
is to assume u s h and B constant during the free expansion 
phase, i.e. up to t — tsr- In order to estimate tsr we con- 
sider two different situations which can represent a typical 
type I/a and a core-collapse supernovae. In the first case the 
supernova explodes in the regular ISM with typical density 
and temperature n\ = 1 cm -3 and To = 10 4 K. On the 
other hand SNRs generated by very massive stars expand 
into medium with higher temperature. This could be either 
the bubble generated by the progenitor's wind or a so called 
super-bubble, a region where an ele vate rate of SN explo- 
sions produces a diluted and hot gas IjHigdon fc Lingenfelter] 
2005). In both cases typical values for the density and tem- 
perature inside the bubble are n\ — 10~ 2 cm -3 and To = 10 6 
K. For both type I/a and core-collapse supernovae we as- 
sume the same value for the explosion energy Esn = 10 51 
erg, and mass ejecta M e j = IAMq. The resulting Sedov 
time is tsr = 470 yr and tsr ~ 2185 yr, respectively. 
The average shock speed during the free expansion phase 
is u s h — 6000 km/s in both cases. 

It is worth stressing that if one would consider a more 
realistic scenario for core-collapse SNe it is important to take 
into account the type of the progenitor and its wind, in order 
to provide a better estimate of the ejected mass and of the 
bubble's dimension and density. These parameters strongly 
affect the value of tsr, which determines how effective is 
the ionization in reducing the maximum energy: indeed for 
tsr 3> 10 3 yr the role of ionization becomes negligible, as 
we see below. 

The last parameter we need is the magnetic field 
strength. We stress that our aim is only to estimate the 
effect of ionization in the scenario where the SNRs indeed 
produce the observed CRs spectrum and the knee is inter- 
preted as the superposition of a rigidity-dependent cutoff of 
different species. For this reason we assume that both sce- 
narios are able to accelerate protons up to the knee energy, 
i.e. E kncc = 3 • 10 15 eV. This condition gives Bi = 160uG 
and 35/j.G for type I/a and core-collapse cases, respectively. 
We note that the chosen values of magnetic field strength are 
consistent with those predicted by the CR-induced magnetic 
filed amplification. 

In Fig. we plot the maximum energy achieved at the 
beginning of the Sedov- Taylor phase by different chemical 
species, from H up to Zn (Zn = 30). The two panels show 
the case of type I/a (upper) and core-collapse SNR (lower) 
previously described. Each panel contains four lines: thin 
solid lines are the maximum energy achieved by ions which 
start the acceleration completely stripped, -Em ax , while the 
remaining lines show -E max computed according to Eq. (|12|l 
for three different locations of the SNR in the Galactic plane: 
at the center of the Galaxy and at 4 and 12 kpc away from 
the Galactic Center. Looking at the upper panel we see that 
the maximum energy achieved by different nuclei in type I/a 
SNRs does not increase linearly with Zn, instead it reaches a 
plateau for Zn ~ 25 at a distance d = 4 kpc and for Zn ~ 15 
at d = 12 kpc. Only for SNRs located in the Galactic bulge 
the reduction of the maximum energy is negligible, at least 
for elements up to Zn = 30. The effect of ionization is much 
less relevant for core collapse SNRs: only those remnants 
located at a distance of 12 kpc show a noticeable reduction 
of J? max . 
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Figure 2. Maximum energy achieved by chemical specie with 
nuclear charge Zjv (ranging from H to Zn) at the beginning of 
the Sedov phase. The upper panel shows the case of type I/a 
SNR, expanding into a medium with density po = 1 cm -3 while 
the lower panel is for core-collapse SNR expanding into a diluted 
bubble (po = 0.01cm -3 ). The thin solid line shows the maximum 
energy achieved by atoms which are fully ionized since the be- 
ginning of the acceleration, while the other curves are computed 
including the photo-ionization due to ISRF for SNR located at 
three different distance d from the Galactic Center, as specified 
in the caption. 





t.ST (yr) 


d (kpc) 


ion ^ J I 


-Emax/^max 


type I/a 


470 





76 (137) 


0.98 (0.52) 






4 


356 (385) 


0.80 (0.39) 






12 


1446 (2932) 


0.48 (0.24) 


core- 


2185 





113 (142) 


0.99 (0.58) 


collapsed 




4 


264 (374) 


0.97 (0.48) 






12 


2860 (2780) 


0.65 (0.35) 



Table 1. Total photo-ionization time and normalized maximum 
energy achieved by iron nuclei for the two idealized cases of type 
I/a and core-collapsed SNRs, which are located at a distance d = 
0, 4 and 12 kpc from the Galactic Center. For comparison also the 
Sedov time is listed. The bare numbers represent the case where 
the shock speed and the magnetic field are both constant during 
the free expansion phase (as in Fig. |2j, while the numbers in 
round brackets are computed assuming and B\ which evolve 
in time according to the model explained in the text (see also 
Fig. 0. 



Figure 3. The same as in Fig. [2] but for a scenario where both 
the shock speed and the magnetic field strength evolve with time, 
as explained in the text. 



The numerical results for iron nuclei corresponding to 
all cases shown in Fig.[2]are summarized in Table[T] Here we 
report the total ionization time and the ratio between the 
energy achieved at tsr and the maximum theoretical energy 
achieved by bare nuclei, fiU lax (notice that the numbers in 
brackets refer to Fig. [3] as explained below). 

The results presented in Fig. [2] are inferred using a sim- 
ple steady state approach. This approach could be too re- 
ductive and one can guess how the effect due to the ion- 
ization process changes when the full evolution history of 
the remnant is taken into account. Needless to say a cor- 
rect calculation of the maximum energy requires not only 
the treatment of the evolution of the remnant but also the 
inclusion of nonlinear effects which consistently describe the 
magnetic field amplification and the back reaction of CRs 
onto the shock dynamics. A fully consistent treatment of 
the problem is beyond the purpose of this work. However we 
would understand whether the evolution and the nonlinear- 
ity can reduce or exacerbate the effect induced by ionization. 
In order to reach this goal we approximate the continuum 
evolution into time steps equal to the ionization times, r p h,k 
and we assume that the stationary approximation is valid for 
each time step. Under this approximation Eq. (|12p can be 
used to include the effect of evolution if we make reasonable 
assumptions on how the shoc k speed and the magnetic field 
change with time. We follow iTruelove fc McKed i|l999l) for 
the description of the time evolution of the remnant. Specif- 
ically we adopt the solution for a remnant characterized by 
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a power-law profile of the ejecta (in the velocity domain) 
with index n = 7, ex panding into an homogene ous medium. 
For this specific case iTruelove fc McKed (|l999l ) provide the 
following expression for the velocity of the forward shock 
(see their Table 7): 

ui(t)/u ch = 0.606 (t/t ch )~ 3/7 (fort < t ST ) (13) 

where u ch = (E SN /M ej f' 2 , t ch = E^J 2 M?/ 6 p" 1/3 and 
tsr ~ 0.732 t c h- We adopt the same values of Esn, M e j and 
pi used in the two cases above. 

For what concern the magnetic field we assume that the 
amplification mechanism is at work, converting a fraction 
of the incoming kinetic energy flux into magnetic energy 
density downstream of the shock. The simplest way to write 
this relation is: 

Bf(t)/(87r) = aBpiUi(t). (14) 

We always assume B2 = 4_Bi. The parameter as hides 
all the complex physics of magnetic amplification and par- 
ticle acceleration. For the sake of completeness we men- 
tion here that in the case of resonant streaming instability 
ob oc ^crVA/ui, where £ cr is the efficiency in CRs and va 
is the Alfven velocity, while in th e case of resonant ampli- 
fication clb oc £ cr iii/4c (Bell 2004). However such relations 
cannot be applied in a straightforward way since they re- 
quire using a non linear theory. Here we prefer to make the 
simplest assumption, taking os as a constant. We determine 
its value from the same condition previously used to fix the 
value of the constant magnetic field: namely we assume that 
for t = tsr the energy of protons is 3 • 10 15 eV. This con- 
dition gives a B = 3.25 ■ 10" 3 and 6.80 ■ 10" 3 for type I/a 
and core-collapse cases, respectively. It is worth noting that 
in the case of some young SNRs the value of as has been 
estimated from the measurement of both the shock speed 
and the magnetic field strength. However we stress that the 
magnetic field cannot be determined in a completely model 
independent way. A possible method requires the measure- 
ment of the spatial thickness of the x-ray filaments: assuming 
that this thickness is determined by the rapid synchrotron 
losses of radiating electrons we can get the average magnetic 
field (see Tableland the listed references and also Table 1 
from lCaprioli et all (|2009l )). Noticeably the values estimated 
for aB are only a factor 5—10 larger then the values we use 
here. 

Now Eq. (|12p can be used to compute the maximum 
momentum at t — tsr- For each time-step, r p h,k, the values 
of Uk and Bt are computed according to Eqs. (|13[) and (|14[) . 
respectively, evaluated at the beginning of the time-step. In 
Fig. [3] we report the results for the maximum energy for 
the case of type I/a (upper panel) and core-collapse SNRs 
(lower panel). The lines have the same meanings as in Fig [21 
As it is clear from both the panels, the reduction of the 
maximum energy is now more pronounced with respect to 
the stationary case shown in Fig [2] Even in the scenario 
of a core-collapse SNR located in the Galactic bulge, iron 
nuclei are accelerate only up to about one half of the max- 
imum theoretical energy. The numerical values of TjJSJj and 
-Emax/£max f° r i ron are reported in Table [T] in round brack- 
ets. One can see that even in the cases where t££ does not 
change with respect to the stationary case, the maximum 
energy results smaller. 

The reason why ionization is less effective when the evo- 




10 20 50 100 200 500 10002000 
t[yr] 

Figure 4. Behavior of maximum energy of Fe ions as a function of 
the time in the case of core-collapse SNR located in the Galactic 
bulge. Two different scenarios are shown: (thin lines) acceleration 
at a constant rate (B\ = const and u\ = const); (thick lines) 
acceleration including the time evolution of the remnant (B\ = 
B\(t) and u\ = u\(t)). The dashed lines represent the case of 
completely ionized ions, while the solid lines include the effect of 
photo-ionization. 



lution is taken into account is a consequence of the fact 
that acceleration occurs mainly during the first stage of the 
SNR expansion. In fact both the shock speed and the mag- 
netic field strength are larger at smaller times. On the other 
hand the effective ions' charge is small during the very ini- 
tial phase of the expansion, hence the acceleration rate is 
smaller than its maximum possible value. In order to clar- 
ify this point in Fig. [4] we compare the energy achieved by 
iron nuclei as a function of the time, for the case of core 
collapse SNR located in the Galactic bulge, with and with- 
out the evolution. The dashed lines shown the case of com- 
pletely ionized nuclei, while the solid lines take into account 
the ionization process. When the acceleration rate is con- 
stant (thin lines) the bulk of the acceleration occurs close 
to tsr, when the ions are almost completely ionized, hence 
i?max and i?^, are practically the same. On the other hand 
when the evolution is taken into account (thick lines) the 
acceleration mainly occurs during the first 200 yr when the 
ionization is still not complete and for t = tsr we have 
-Emax/£max = 0.58. Besides the numerical value, this exer- 
cise shows that the inclusion of the time evolution and non 
linear effects can enhance the reduction of the maximum 
energy due to the ionization process. 

Some comments are in order. The most relevant conse- 
quence of the ionization mechanism concerns the shape of 
the CR spectrum in the knee region. As we have already dis- 
cussed, the relation -E ma x,jv oc Zn is needed in order to fit 
the data. Even a small deviation of the cutoff energy from 
the direct proportionality can affect the slope. Let us as- 
sume that the measured slope, s = 3.1, indeed results from 
this proportionality relation. We have i5 m ax,Fe = 26E maX)P . 
Now if the maximum energy of iron decreases to 80% of its 
maximum value, the slope changes to ~ 3.33, while for a 
decrease of 50% the slope becomes ~ 3.94. 

According to our results we can say that in the context 
of the SNR paradigm, the primary sources of CRs above the 
knee energy are most probably the core-collapse SNRs with 
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M e j ^> IMq. Type I/a SNRs seem unable to accelerate ions 
up to an energy Zn times the proton energy, due to their 
small Sedov- Taylor age. This consideration is strengthened 
by the fact that the CR flux observed at the Earth is mostly 
due to the SNRs located in the solar neighborhood, rather 
than those located in the bulge. In fact the escaping length 
from the Galaxy is determined by the thickness of the Galac- 
tic halo which is ~ 3 — 5 kpc, while the Galactic bulge is at 
8 kpc from us. 

A second comment concerns the acceleration of ele- 
ments beyond the iron group. Even if the contribution of 
such ultra-heavy elements to the CR spectrum is totally 
negligible at low energies, at higher energies it can hardly 
be measured. In principle the contribution of ultra-heavy el- 
ements could be significant in the 100 PeV regime, which is 
the energy region where the transition between Galactic and 
extragalactic CRs occurs. Indeed several authors pointed out 
that a new component is needed to fit this transition region, 
beyond the elements up to iron accelerated in "stan dard" 
SNRs (see e.g. the di s cussio n in ICaprioli et al.1 {2010)). As 
inferred by iHorandell (|2003h stable elements heavier than 
iron can significantly contribute to the CR spectrum in the 
100 PeV regime if one assume that their maximum energy 
scales like Zn- This assumption is especially appealing also 
for a second reason: in principle it could explain the pres - 
ence of the second knee in the CR spectrum (Horandcl 2003); 
in fact the ratio E 2 „d kn< , e /Eknee — 90 is very close to the 
nuclear charge of the last stable nucleus, uranium, which 
should have -Emax.fj = 92_E max , p = 414 PeV. On the other 
hand the ionization mechanism discussed here provides a 
strong constraint on the role of ultra-heavy elements. If the 
acceleration of elements heavier than Fe occurs at SNRs like 
those considered above, it is easy to show that they can- 
not achieve energies much larger than the Fe itself, because 
the total ionization time increases rapidly with the nuclear 
charge. A contribution from ultra-heavy elements is maybe 
possible only if they are accelerated in very massive SNRs, 
with M e j S> IMq. A second possibility that would be in- 
teresting to investigate is the acceleration during the very 
initial stage of the remnant, when the expansion occurs into 
the progenitor's wind. In this case the high density of the 
wind could strongly reduce the ionization time thanks to the 
Coulomb collisions and nuclei could be totally stripped in a 
short time. 



4.1 Spectrum of ions 

In order to compute the distribution function of a single 
chemical specie we consider separately each ionization state. 
Let /jv(£iP) be the distribution function of the specie N 
with effective charge Z. The equation that describe the dif- 
fusive transport of ions in one dimension, including the ion- 
ization process reads: 

df N d 2 f N idu df N f N z 
u ^x~ = D{p) ^ + 3 d7 p ^p~ ~ 7W) Qn{x ' p) ' (15) 

where u{x) is the fluid velocity and D(jp) is the diffusion 
coefficient, which is a function of Z and it is assumed con- 
stant in space. Tjf(p) is the ionization time for the losses 
of one electron, namely for the process N z — > N z+1 + e~ . 
We neglect multiple ionization processes. In order to solve 
Eq. (|T5|) we need to specify the injection term Q N (x,p). We 
assume that ions with the lower degree of ionization, Zo, are 
injected only at the shock position (x = 0) and at a fixed 



momentum pi n j . Hence the source term for is 
Ql?{x,p)=K N 5{x)8(p-p- mi ). 



(16) 



The normalization constant K N ° is determined from the to- 
tal number of particles injected per unit time. For the sub- 
sequent ionization states, Z > Zo, the injection comes from 
the ionization of atoms with charge equal to Z — 1, i.e.: 



x,p) 



/jv(^,p)/tjv(p) • 



(17) 



Equation (|15[) can be solved separately in the upstream and 
downstream regions, where the term du/dx vanishes and 
the equation reduces to an ordinary differential equation 
of second order. We label the quantities in the upstream 
(downstream) with a subscript 1(2) and we adopt the con- 
vention x < (>)0 in the upstream (downstream). We fix 
the boundary conditions at the shock position such that 
f(x = 0,p) = fo(p), where /o(p) has to be determined. The 
boundary condition at infinity is df/dx = for x — > ±oo. 
For the sake of simplicity we now drop the label character- 
izing the ion specie, such that /jf = /, and we write the 
solution in a compact form as follows: 



fi(x,p) = fo(p)e^T* ± 



Qi(x ,p) e ±Xj±(x-x') j^i 
iVl + A t 

Qi(x ,p) gT^Ca-gQ^/ 



o Ui VI + Aj 



/ Qi{ x \p) -Xi-I 



Wl + Ai 



dx . (18) 



4 SOLUTION OF THE TRANSPORT 
EQUATION 

In this section we present the solution of the stationary 
transport equation for ions and electron accelerated at 
shocks when the ionization process is taken into account. 
We first get the solution for the distribution function of ions 
which can be used, in turn, as the source term responsible 
for the injection of electrons. We will consider only the test- 
particle solution, neglecting all the non linear effects coming 
from the back reaction of the accelerated particles onto the 
shock dynamics. 



The upper and lower signs refer to the downstream (i — 2) 
and upstream (i = 1) solutions, respectively. The dimension- 
less function A(p) is the ratio between the diffusion time and 



the ionization time, i.e. A;(f>) 



u' 2 r(p) ' 



\i+ and Xi- are 



the inverse of the propagation lengths along the counter- 
streaming and the streaming directions with respect to the 
fluid motion, respectively, i.e.: 



2Di(p) 



*h/l + Ai(p)±l] 



(19) 



The quantity A takes into account the diffusion, the ad- 
vection and the ionization processes. One can see that in 
the limit where the ionization is negligible, namely when 



The role of ionization in the shock acceleration theory 9 



t — > oo, the streaming propagation length diverges, while 
the counter-streaming one reduces to D/u. 

The three integral terms which contribute to the distri- 
bution function in Eq. (I18|l have a clear physical interpre- 
tation. Let us consider the solution in the upstream: in this 
case the first integral represents the contribution due to par- 
ticles advected with the fluid from — oo, the second accounts 
for the particles diffusing in the counter-streaming direction, 
while the last integral takes into account the flux escaping 
across the shock surface. For the downstream solution the 
meanings of the first end second integrals are reversed, while 
the third one keeps the same meaning. This interpretation 
becomes clear if one take the limit r — > oo. 

From Eq. (|18jl we can easily recover the standard solu- 
tion where the ionization is not taken into account, assuming 
that the injection occurs only at the shock position and us- 
ing Aj = 0: in this case the integral terms vanish, and we 
find fi(x,p) = / (p)exp[-ttia;/£)i] and fn(x,p) = /o(p). 

In order to get the complete solution we need to com- 
pute the boundary distribution function at the shock posi- 
tion, /o(p). We follow a standard technique which co nsists 
in int egrating Eq. (|15[) around the shock discontinuity (Blasi 
2002). We got the following differential equation for fo(j>)' 

=( D >f)„ + (^)„-' <*> 

where + and 0~ indicate the positions immediately down- 
stream and upstream of the shock. The quantities Ddf /dx 
can be easily computed deriving / with rspect to x from 
Eq. JIB]). We get: 

(^f ) 0+ = f °Y ( 1 + ^^jf V(21) 



and 




After a little algebra, the solution of Eq. (|20[) can be ex- 
pressed in the following form: 

Mp) = s P -° [ P ^/^exp{- S r^h( P ")\ , (23) 

where the usual definition of the power law index is s — 
3mi/(ui —112). The injection term G is the sum of the up- 
stream and downstream contributions, 

/0 poo 
Q ie Xl - x 'dx'+ / Q 2 e- X2 + X 'dx' (24) 
00 J 

while the function h in the exponential is: 

hip) = \Wl + A 1 (p) - 1) + 1(0 + A 2 (p) - 1) . (25) 

The solution in Eq. (|23|) shows the typical power law behav- 
ior oc p~ s in the region of momentum where neither injection 
nor ionization are important, in fact when A(p) <C 1 then 
h(p) vanishes. On the other hand when the ionization length 
becomes comparable with the diffusion length, i.e. A(p) ~ 1, 
the function h(p) produces an exponential cutoff in the dis- 
tribution function. 

Now that we got the formal solution we can show some 
results. Let us start considering the acceleration of He. We 
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Figure 5. Distribution function of Helium at the shock position, 
/hc o(p)- The thin solid lines represent the distributions of He+ 
and He++, while the thick line is the sum of both the contri- 
butions. The dotted lines represent the distribution of electrons 
(multiplied by 10 4 ) ejected in the process He + — > He+" 1 " + e — , 
and subsequently accelerated. 



assume that the initial ionization state at the moment of in- 
jection is He + . The distribution function / Hc + is computed 
with the procedure described above, using the injection term 
in Eq. (|16[1 . We fix p- ln j — 10~ 3 m p c and the shock velocity 
equal to 3000 km/s. The photo-ionization time can be an 
arbitrary function of p. We use the approximate expression 
given in Eq. ©. Eqs. (|18|l and (I23|l can be integrated nu- 
merically. Once / He + is known, we can compute /h c ++ using 
the injection term in Eq. ([T7)l . i.e. with Q Hc ++ — /h c + / r Hc+ • 
The resulting distributions are plotted in Fig. [S] with thin 
solid lines. We note that / Ho + x P 4 up to p ~ 10m p c; above 
such value the distribution of He + drops exponentially while 
the distribution of He ++ starts rising. It is worth noting 
that the total distribution (plotted with thick solid line) is 
slightly steeper than p~ 4 in the transition region. This effects 
is due the different probability which He + and He ++ have to 
return at the shock once they are downstream. In fact an ion 
with charge Z located at the position x in the downstream 
have a probability to return at the shock oc e~ U2X ^ D( - z ^; after 
the ionization the charge becomes Z+l and the return prob- 
ability is reduced by a factor e~^ z+1 ^ z . This translates into 
a net loss of particles from downstream. This effect of losses 
is more pronounced when considering heavier elements. In 
Fig. [6] we show the case of carbon atoms which start the 
acceleration as singly ionized. The distribution functions of 
all ionized state are plotted with thin solid lines, from C + 
up to C 6+ . In the energy region where ionization occurs, 
1.5 < 7 < 4.5, the slope of the total distribution function is 
s ~ 4.1. 

4.2 Spectrum of electrons 

Now that we know how to compute the distribution function 
of ions, the solution for the electron distribution is straight- 
forward. We have to solve the same transport equation (|15|l 
writing down the right expression for the electron injection, 
Q e , and dropping the term due to ionization. In the present 
calculation we neglect energy losses. 

Let us consider a single specie N. From each photo- 
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Figure 6. Distribution function of Carbon in all the ionization 
state from C + up to C 6+ . The thick solid line is the sum of all 
contributions. The spectra of electrons ejected from each ioniza- 
tion process are plotted with thin dashed lines, while the thick 
dashed line represent the total sum (xlO 4 ). 



ionization process of the type N z — > N z+1 + e , we have 
the following contribution to the electron injection: 

Qf z {x,p) = r d p'c^4r6^{p-^p') 

J p T N(P ) 



= 6 



\fN,l( X 'P/£ N ) 



(26) 



The 5-Dirac function in momentum is present because elec- 
trons are ejected with the same Lorentz factor of the parent 
ions, as we discuss in Ej2j Hence the electron momentum is 
Pe = PNirLe/rriN = P n£,n- Following the procedure used in 
the previous section, we can write the distribution function 
for electrons at the shock position as follows: 

fP '^^. (27) 

v u i 

Pin] 



rN.Z 1 \ -S/-S-3 
fe,0 (P) = S P CiV 



Equation (|27|l is similar to Eq. (|23|l with the exception of 
the multiplying factor and for the absence of the ex- 

ponential cutoff. We note that the contribution to electron 
injection from the downstream is negligible with respect to 
the upstream. This occurs because electrons produced in the 
downstream have a negligible probability to come back to 
the shock. The return probability for one electron stripped 
downstream at a position a; is a factor e ~ m «/ z «"ic sma vj er 
that the return probability of the parent ion. Hence in the 
injection term we consider only the contribution from up- 
stream, which reads: 



''(jp) = 



' z (x',p)dx' . 



(28) 



In Figs. [S] and [51 besides the distribution of ions, we plot 
also the distributions of the electrons due to each single ion- 
ization step, multiplied by 10 4 (dashed lines). Comparing 
Eqs. (|23p and (|27[1 we see that the ratio between electron 
and ion distribution functions is: 



KeN = / e A Z (P)//iV,o(p) = £nIn,z(p) ■ 



(29) 



The function In,z{ P ) represents the ratio between the in- 
jection integrals in Eqs. I|23[) and (|27[) . It can be computed 



numerically and in the asymptotic limit p>p* (where p* 
is such that A(p*) = 1) it is I N ,z ~ Z/(2Z - 1). This result 
is not surprising: it reflects the fact that for each ionization 
event electrons can be injected only if they are released up- 
stream hence their number is roughly 1/2 of the injected 
ions. 

Now we are interested in the total number of accelerated 
electrons that can be produced via the ionization process. In 
the literature the number of electrons is usually compared 
with that of accelerated protons: DSA operates in the same 
way for both kinds of particles, hence a proportionality rela- 
tion between their distribution functions is usually assumed, 
i.e. f e (p) = K ep f p (p). From the observational point of view, 
the value of K ep in the source strongly depends on the as- 
sumption for the magnetic field strength in the region where 
electrons radiate and, in the context of the DSA theory, it 
can be determined for those SNRs where both non thermal 
X-ray and TeV radiation are observed. Two possible scenar- 
ios have been proposed. In the first one electrons produce 
both the X-ray and the TeV components via the synchrotron 
emission and the inverse Compton effect, respectively; this 
scenario requires a downstream magnetic field around 20/iG, 
and K ep ~ 1CT 2 - 10" 3 . The second scenario assumes that 
the number of accelerated protons is large enough to ex- 
plain the TeV emission as due to the decay of neutral pions 
produced in hadronic collisions. In this case the DSA re- 
quires a magnetic field strength of few hundreds fiG and 
K ep ~ 10 -4 — 10 -5 . Such a large magnetic field is consis- 
tently predicted by the theory as a result of the magnetic 
amplification mechanisms which operate when a strong CR 
current is present. In Table [2] we report some examples of 
young SNRs where the nonlinear theory of acceleration has 
been applied to fit the multi- wavelength spectrum. The forth 
and fifth columns report the estimated valued of the down- 
stream magnetic field strength and the corresponding elec- 
tron/proton ratio. 

In order to give an approximated estimate for K ep we 
need to multiply Eq. (|29|) by the total number of ejected 
electrons, i.e. (Zn — Zn,o), and sum over all atomic species 
present in the accelerator, i.e. K ep ~ ^2 N Kn p (Zn — 
Zjv i6 ff) KeN ■ Kn p are the abundances of ions measured at 
the source in the range of energy where the ionization oc- 
curs. Even if the values of Kn p are widely unknown, wc 
can estimate them using the abundances measured at Earth 
and adding a correction factor to compensate for propaga- 
tion effects in the Galaxy, namely the fact that particles with 
different Z diffuse in a different way. The diffusion time in 
the Galaxy is usu ally as s umed to be t^h oc (p/Z)~ s , with 
5 » 0.3 - 0.6 fsee iBlasil (|2007T ) for a review on recent CR 
experiments). If Kn p ,o is the ion/proton ratio measured at 
Earth, than the same quantity measured at the source is 



K 



Np 



K 



Np,0^N 



Hence the final expression for the elec- 



tron/proton ratio at the source is: 



Kep = K 
N 



n p ,o Z N (Zn ■ 



j Zn / m e 
2Zn — 1 \rriN 



.(30) 



In Fig. [7] we report the value of K ep computed according 
to Eq. (|30|) as a function of the spectral slope s and for 
5 — 0.3 and 0.6. For Kn p ,o we use the abundances of nuclei 
in th e CR spectrum measured at 1 TeV ( Wiebel-Sooth et al.l 
1998). Moreover we determine the initial charge of different 
species, Zn,o, adopting the thermal equilibrium in a plasma 
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SNR 


Age (yr) 


d(kpc) 


B 2 (^G) 


K ep /10~ 4 


rcf. 


Kepler 


400 


3.4-7.0 


440-500 


0.7-2.8 


(a) 


Tycho 


430 


3.1-4.5 


368-420 


4.2-15 


(&) 


SN 1006 


1000 


1.8 


90 


1.3 


(c) 






2.2 


120 


1.1 


(d) 


G347.3-0.5 


1600(?) 


1.0 


100 


0.6 


(e) 






1.0 


130 


1.0 


(/) 



Table 2. Estimated value for the magnetic field and K ep for some 
young supernova remnants. Results are from: ( a 'Berezhko et al. 
(2006), ^Volket al. (2008), ( c >Morlinoet al. (2010), ^Berezhko 
et al. (2009), ( e )Morlinoet al. (2009), ^Berezhko & Volk (2006). 
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Figure 7. Value of K ep computed from Eq. (1 30 P> as a function of 
the slope s and for 8 = 0.6 (solid line) and 8 = 0.3 (dashed line). 

with a temperature T ~ 10 5 . Remarkably using the slope 
predicted by linear theory for strong shock, s = 4, we have 
A'ep ~ 10~ 4 : this number gives the right order of magnitude 
required to explain the values of K ep reported in Table [2] 
We stress here that in order to provide a better estimate of 
K ep for a specific SNR we should know the local chemical 
composition and have a realistic model for the injection of 
heavy elements into the DSA, which can provide the right 
value of Zn.o- 

Equation (|30p has been obtained using the test particle 
approach. A better estimate requires the use of nonlinear 
theory which is, indeed, much more complicated to develop. 
Nevertheless it is easy to realize that the main effect that 
nonlinearity can produce on K ep is due to the different pre- 
diction of the slope s. In fact, while the test particle ap- 
proach gives the universal slope s — 4, the nonlinear theory 
predict concave spectra with s = s(p). The basic formulation 
of the nonlinear theory usually leads to s > 4 for momenta 
lower that ~ WZ -1 G eV/c and s < 4 for larger momenta 
l|Amato fc Blasill2005l ). As discussed in i|2] ionization typ- 
ically occurs for 7 > 10, where the slope is s < 4. Hence 
we can conclude that the nonlinear effects tend to increase 
the value of K ep . For instance, from Fig. [7] we see that for 
s — 3.7, K ep rises up to ~ 10~ 3 . 

It is worth noting that nonlinear effects can also result 
in the opposite situation with K ep < 10~ 4 , as we will show 
below. A key ingredient of the DSA is the speed of the mag- 
netic turbulence which is responsible for particles scattering. 
The speed of the scattering centers contributes to the deter- 



mination of the effective compression ratio felt by particles 
which, in turns, determines the spectral slope. The typical 
speed of the turbulence is of the order of Alfven speed which, 
when is computed in the background magnetic field, is of 
the order of few tens km/s and it is negligible with respect 
to the shock speed. On the other hand if one assume that 
turbulence moves with the Alfven s peed computed in th e 
amplified magnetic field, as shown bv lCaprioli et all l|20ld ). 
it cannot be neglected and the resulting slope can be > 4. If 
this happens the value of K ep drops and the number of elec- 
trons produced via ionization is maybe insufficient in order 
to explain the synchrotron emission. 

A final comment concerns the mechanism of electrons 
injection in core-collapse SNRs. The value of K ep shown in 
Fig.[7]has been computed using the ionization level of chem- 
ical elements in a plasma with T ~ 10 s . We recall that the 
type I/a SNRs expand in the warm ISM, where the typical 
temperature is around 10 4 K hence we expect the injection 
of electrons to be relevant. On the other hand core-collapse 
SNRs expan ds into a diluted bubble w hose temperature 
reach 10 6 K l|Higdon fc Linge nfclter 200l[), hence the degree 
of ionization of atoms is higher and the num ber of electrons 
availa ble for the injection is smaller (see e.g. IPorquet et al.l 
i|200ll )). Using T = 10 6 we estimate that K ep is a factor ~ 4 
smaller than the values shown in Fig. [7] On the other hand in 
those bubbles the met allicity is estimated to be grea ter than 
the ISM mean value (|Higdon fc Li ngenfclter 2005|). Hence, 
even for core-collapse SNRs the electron injection through 
ionization could play a relevant role. 



5 DISCUSSION AND CONCLUSIONS 

In this work we include for the first time the ionization pro- 
cess of heavy ions in the theory of DSA at SNR shocks. 
We showed that, for the typical environments where SNRs 
propagate, the photo-ionization due to the ISRF dominates 
on the collisional ionization and produces two important ef- 
fects: 1) the reduction of the maximum energy achieved by 
ions and 2) the production of a relativistic population of 
electrons. 

(1) The first effect is especially important for the inter- 
pretation of the knee structure in the all particle spectrum of 
CRs. The change of slope on the two sides of the knee is gen- 
erally interpreted as due to the superposition of spectra of 
chemicals with different nuclear charges combined with their 
abundances and convolved with the rigidity dependent dif- 
fusion in the Galaxy. An important requirement for a good 
fit to the knee is that the cutoff energy of each specie, from 
hydrogen up to iron, should be proportional to the nuclear 
charge, Zn- 

It is generally assumed that DSA can accelerate ions 
up to a maximum energy proportional to Zn- This assump- 
tion is valid only if the ionization time needed to strip all 
electrons from atomic orbitals, Tjjjjf, is much shorter than 
the acceleration time. We assume that the maximum pos- 
sible acceleration time is equal to the Sedov- Taylor time, 
tsT, corresponding to the end of the free expansion phase. 
Using a simple steady-state approach we showed that, es- 
pecially for very massive ions, r*^ can be comparable to 
or even larger than tsr, depending on the type of the rem- 
nant and on its location in the Galaxy. In fact the ISRF 



12 G. Morlmo 



responsible for the photo-ionization decreases for increasing 
distance from the Galactic Center. We analyze two differ- 
ent cases which are representative of a type I/a and a core- 
collapse SNRs. Our main results are the following. SNRs of 
type I/a are generally unable to accelerate ions according to 
£max,jv oc Zn because tsr does not exceed the typical value 
of 500 yr. The only exception are the remnants located in 
the Galactic bulge where the photon field is strong enough 
to reduce the photo-ionization time to a value much shorter 
than tsr- The situation for core-collapse SNRs is different 
because their typical Sedov time is Z 2000 yr hence the ion- 
ization time for chemicals up to iron can be neglected and 
the maximum energy is indeed proportional to the nuclear 
charge, with the exception of those remnant located very far 
away from the Galactic Center where the ISRF is very low. 

Previous conclusions are based on the assumption that 
during the free expansion phase both the magnetic field and 
the shock velocity remain constant. This is indeed a poor 
approximation. In fact during the free expansion phase the 
shock speed can vary by a factor of few (depending on the 
velocity profile of the ejecta and on the density profile of 
the external medium). Now for DSA the acceleration rate 
is proportional to u 2 sh SB(u B h), where SB(u s h) is the turbu- 
lent magnetic field generated by the CR induced instabilities 
(resonant or non-resonant), which is an increasing function 
of the shock velocity. Hence the acceleration rate can signif- 
icantly change during the free expansion phase even if u B h 
varies only by a factor of few. Always in the framework of 
steady-state approximation, we develop a simple toy model 
able to include the magnetic field amplification and the time 
evolution of the remnant. Using this toy model we showed 
that the magnetic amplification and time evolution can sig- 
nificantly reduce the acceleration time with respect to tsr 
and, as a consequence, the maximum energy achieved by 
heavy ions becomes smaller than the previous estimate: the 
relation £ mlD . ; jv oc Zn can be achieved up to iron nuclei only 
by SNRs with a Sedov time ^> 10 3 yr which means that the 
mass ejecta should be 3> I Mq. 

Clearly the same mechanism also applies to ions heavier 
then Fe. Their maximum energy cannot be much larger than 
that achieved by Fe itself unless the acceleration occurs in 
very massive SNRs. Hence the ionization mechanism put 
severe limitation to the possibility that ultra heavier ions 
can contribute to the CR spectrum in the transition region 
between Galactic and extragalactic component. 

(2) The second effect concerning the pr oduction of rela - 
tivistic electrons was already put forward in iMorrmo^ OQcn. 
The ionization mechanism provides a source for the injec- 
tion of relativistic electrons into the DSA mechanism. We 
investigate the possibility that those electrons can be re- 
sponsible for the synchrotron radiation observed from young 
SNRs. In order to estimate the total number of injected elec- 
trons we use a semi-analytical technique to solve the steady- 
state transport equation which describe the electron distri- 
bution function in the shock region. Summing the contribu- 
tions coming from the ionization of all chemicals we showed 
that the ratio between accelerated electrons and protons is 
K ep ~ 10~ 4 . This value is especially appealing because cor- 
responds to the right order of magnitude required in order to 
explain the synchrotron emission if the magnetic field is am- 
plified up to few hundreds ^iG. Interestingly such magnetic 
field amplification is consistently predicted by non linear ac- 



celeration theory as a consequence of the back reaction of 
accelerated particles onto the shock dynamics. 

Several effects can modify our prediction for the num- 
ber of electrons. The estimate K ep ~ 10~ 4 is obtained using 
test particle approximation and adopting the abundances of 
CR spectrum observed at the Earth corrected for the effect 
of propagation in the Galaxy. Even if we did not develop a 
fully nonlinear approach, we showed how nonlinear effects 
can substantially modify the electron/proton ratio. Indeed 
the value of K ep strongly depends on the spectral slope as 
one can see from Eq. ()30[) . Nonlinear DSA usually predicts 
steeper spectra than the test particle result, s — 4, and this 
translates into an enhancement of K ep . On the other hand 
also the opposite situation can be realized. Indeed some au- 
thors pointed out that, if the magnetic field amplification 
occurs, the velocity of the magnetic turbulence with respect 
to the plasma can be much larger than the Alfven speed cor- 
responding to the background magnetic field . If this occurs 
the spectral slope of the accelerated particles be apprecia- 
bly softer than 4 and the electron/proton ratio would drop 
to a value which is insufficient to explain the observed syn- 
chrotron radiation. 

In order to get a correct determination of K ep for a sin- 
gle SNR, other effects should be taken into account: a) the 
chemical composition of the environment where the SNR ex- 
pands, b) the initial ionization state of each chemicals and 
c) a realistic model for the ions injection. A further com- 
plication directly related to the previous points is the fact 
that many heavy elements are condensed in solid dust grains 
both in the ISM as well as in the stellar wind. It has been 
pointed out that the injection of refractory elements into the 
DSA can be dominated by the sputtering of these dust grains 
rather tha n by the injection of single at oms from the ther- 
mal bath (lEllison. Drurv fe Meverl 0-9971 ). In this picture a 
correct computation of the number of injected ions and their 
initial ionization state is very challenging. 

As a finally remark we want to stress that the elec- 
tron/proton ratio in a single source, what we call K ep , is dif- 
ferent from the same ratio measure d from the C R spectrum 
at Earth, which is 5- 10" 3 at 1 TeV (|Brasill2007T ). Sometimes 
these quantities are assumed to be the same. Conversely 
they could be different because the latter is the sum of the 
contribution coming from all sources integrated during the 
source age, and also reflect transport to Earth and losses in 
transport. Especially relevant is the fact that other sources 
like pulsar wind nebulae can significantly contribute to the 
electron flux measured at Earth, but not to the proton flux. 
Moreover the value of K ep in a single source can vary during 
the source age. Different mechanisms of electron injection, 
other than the ionization, can be relevant in different phases 
of the remnant. 
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